CROSS-SECTIONS FOR SQUARK AND GLUINO 
PRODUCTION AT HADRON COLLIDERS 



R. HOPKER 

Deutsches Elektronen- Synchrotron DESY, D-22603 Hamburg, Germany 

W. BEENAKKER 
Instituut-Lorentz, University of Leiden, The Netherlands 

We present the cross-sections for the hadroproduction of squarks and gluinos 
in next-to-leading order of supersymmetric QCD. The four possible final states 
squark-antisquark, squark-squark, gluino-gluino and squark-gluino are analysed 
for the hadron colliders Tevatron and LHC. The dependence of the cross-sections 
on the renormalization and factorization scale is reduced significantly. The shape 
of the transverse-momentum and rapidity distributions remains nearly unchanged 
when the next-to-leading order SUSY-QCD contributions are included. The size 
of the corrections at the central scale, given by the average mass of the produced 
particles, varies between -|-5% and +90%. 



1 Introduction 

The colored squarks {qL,QR) and gluinos (g), the supersymmetric partners of 
the quarks (q) and gluons (5), can be searched for most efficiently at high- 
energy hadron colliders. As R-parity is conserved in the QCD sector of the 
Minimal Supersymmetric Standard Model (MSSM), these particles are always 
produced in pairs. At the moment they can be searched for at the Fennilab 
Tevatron, a pp collider with a centre-of-mass energy of l.STeV [Ref.Bo]. In 
the future the CERN Large Hadron Collider (LHC), the pp collider with an 
envisaged centreyof-mass energy of 14TeV, will allow to cover mass values up 
to 1-2 TeV [ReiBi. 

So far most of the experimental analyses have been based on the lowest- 
order (LO) production cross-sections. For obtaining adequate theoretical pre- 
dictions the LO cross-sections are in general not sufficient. The most important 
arguments in favor of an analysis that takes into account the next-to-leading- 
order (NLO) SUSY-QCD corrections arc: 

• The LO cross-sections have a strong dependence on the a priori unknown 
renormalization scale Qji and the factorization scale Qp. Consequently, 
the theoretical predictions have in general an uncertainty that is almost 
as large as the cross-section itself. By implementing the NLO corrections 
a substantial reduction of the scale dependence is expected. 
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• From experience with similar processes (e.g. hadroproduction of top 
quarks), the NLO QCD corrections are expected to be sizeable. 

• An enhancement of the cross-section would lead to a higher value for the 
lower mass bounds for squarks and gluinos. 

• In case of discovery of squarks and gluinos, a precise knowledge of the 
total cross-sections would be mandatory for the determination of the 
masses of the particles. 

Here we report on the calculations of NLO SUSY-QCD corrections to the 
production of squarks and gluinos in pp/pp collisions, based on the studies 
presented in Ref.u. For a more comprehensive report on this topic we refer to 
Ref.y. 



2 Technical set-up 

We consider the following hadronic production processes: 

pp/pp qq, qq, gg, qg {q ^ i), (1) 

where the chiralities and flavors of the squarks (e.g. ul, dn) as well as the 
charge-conjugate final states (e.g. qq) are implicitly summed over. They are 
exemplified in Fig.|l| for squark-gluino production. 




Figure 1: Generic Feynman diagrams for squark-gluino production in pp/pp collisions. 

In analogy to the experimental analysis, we exclude top-squarks from 
the final state and take all produced squarks to be mass degenerate. A 
study of the production of pairs of top-squarks is in progress, including the 
(model-dependent) mixing effects in the squark sector. At the partonic level 
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many subprocesses contribute at LO and NLO, corresponding to different fla- 
vors/cfiiralities of the squarks and different initial-state partons. The initial- 
state partons are made up of the massless gluons and the five light quark flavors 
(n f = 5) , which are considered to be massless as well. 

The NLO SUSY-QCD corrections include the virtual corrections (consist- 
ing of self-energy corrections, vertex corrections Fig.H(a), and box diagrams 
(b)), real-gluon radiation (c), and the radiation of a massless quark (d). 




(a) (b) (c) (d) 



Figure 2: Generic virtual corrections [(a): vertex contributions, (b); box contributions] and 
real radiation [(c); gluon radiation, (d): quark radiation]. 

For the particles inside the loops we use the complete supersymmetric 
QCD spectrum, i.e., gluons, gluinos, all quarks, and all squarks. We have 
excluded the top-squarks from the final states. In order to have a consistent 
NLO calculation, however, we have to take the top-squarks into account inside 
loops. For the sake of simplicity we take them mass-degenerate with the other 
squarks. For the top quark we use rrit = 175 GeV. Consequently the final 
results will depend on two free parameters: the squark mass vriq and the giuino 
mass Trig. 

The divergences appearing in the NLO corrections are regularized by per- 
forming the calculations in n = 4 — 2e dimensions. These divergences consist 
of ultraviolet (UV), infra-red (IR), and coUinear divergences, and show up as 
poles of the form {i = 1,2). For the treatment of the 75 Dirac matrix, en- 
tering through the quark-squark-gluino Yukawa couplings, we use the 'naive' 
scheme. This involves a 75 that anticommutes with the other gamma matrices. 
This is a legitimate scheme at the one-loop level for anomaly-free theories. 

The UV divergences can be removed by renormalizing the coupling con- 
stants and the masses of the massive particles. In the case of the mass renor- 
malization we have used the on-shell scheme and for the renormalization of the 
QCD coupling constant the modified JVIinimal Subtraction [MS ) scheme. To 
preserve the supersymmetric relation between the gauge coupling (e.g. in the 
quark-quark-gluon vertex) and the Yukawa coupling (e.g. in the quark-squark- 
gluino vertex) we have to add a finite renormalization to the Yukawa coupling 
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. The heavy particles (top, squarks and giuinos) are decoupled from the 
running of the strong coupling. 

After carrying out this renormalization program, the cross-sections are UV 
finite. Nevertheless there are still divergences left over. The IR divergences 
cancel in the sum of virtual corrections and soft-gluon radiation. In order to 
separate soft from hard radiation a cut-off A is introduced in the invariant mass 
corresponding to the radiated gluon and one of the produced massive particles. 
If soft and hard contributions are added up, any A dependence disappears from 
the cross-sections for A ^ 0. The remaining coUinear singularities, finally, can 
be absorbed into the renormalization of the parton densities, carried out in the 
MS mass-factorization scheme. 

If the giuinos are lighter than the squarks, giuinos can also be decay prod- 
ucts of on-shell squarks, q ^ gq [see e.g. Fig.||(d)]. However, these situations 
are already accounted for by the LO cross-section (e.g. qq production) and the 
subsequent decay. In order to avoid double counting, we restrict ourselves to 
irreducible final states in which giuinos do not evolve from on-shell squarks. 
For the wedge rrig > rriq we disregard, in the same sense as above, the decay 
of the giuinos to squarks. 

3 Results 

The hadronic cross-sections are obtained by the convolution of the parton 
densities with the corresponding partonic cross-sections. When discussing LO 
and NLO results, we calculate all quantities [as{Q^), the parton densities, and 
the partonic cross-sections] in LO and NLO, respectively. We exemplify our 
results by gluino-pair production at the Tevatron and squark-gluino production 
at the LHC. The qualitative behavior of the other processes is quite similar. 

As shown in Fig.|[ we find that the theoretical predictions for the pro- 
duction cross-sections are nicely stabilized by taking into account the NLO 
SUSY-QCD corrections. In all processes, for both the Tevatron and the LHC, 
it is observed that the dependence on the renormalization/factorization scale 
Q {Q = Qr = Qf) is quite steep and monotonic in LO, whereas the Q de- 
pendence is reduced significantly in NLO. Even a broad maximum develops 
at scales near one third of the average mass of the final-state particles. The 
variation of the cross-sections as a result of different NLO parametrizations of 
the parton densities is < 10% at the Tevatron and < 13% at the LHC, where 
the gluon densities are more important. 

From now on we use GRV94 parton densities and conservatively take as 
default scale Q the average mass of the produced particles. The ii'-factors, 
K = UNLO / cfLO , depend strongly on the process. This is exemplified in 
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Figure 3: Scale dependence of the total hadronic cross-sections for (a): gluino-pair production 
at the Tevatron (rrig = 280GeV, rrig = 200GeV, VS = l.STeV), and (b): squark-gluino 
production at the LHC (m^ = 600GeV, nig = SOOGeV, \/S = 14TeV). Parton densities: 
GRV94 (solid), CTEQ3 (dashed), and MRS A' (dotted). 



Fig.^(a) for the Tevatron and in (b) for the LHC. For both coUider types the 
NLO corrections are positive and large (up to +90%) for the dominant pro- 
duction cross-sections, involving at least one gluino in the final state. The 
corresponding iiT-factors also exhibit a sizeable dependence on the squark and 
gluino masses (especially the rriq dependence of Kgg is large). The NLO cor- 
rections for squark final states are moderate (< +30%). In case of discovery 
of squarks and gluinos, the inclusion of the NLO SUSY-QCD corrections is 
needed for an accurate determination of the cross-sections and the masses. 

Apart from the total cross-sections, also distributions with respect to the 
rapidity (y) and transverse momentum (pt) of one of the outgoing massive 
particles can be studied. The differential cross-sections are shown for gluino- 
pairs at the Tevatron in Fig.^and for squark-gluino at the LHC in Fig.^. They 
are normalized to the total cross-section in the corresponding order. The K- 
factors for these distributions are independent of y for all practical purposes 
and hardly depend on pt, where the NLO corrections render the distributions 
somewhat softer. Consequently, multiplication of the LO distributions with 
the above-defined if-factors for the total cross-sections approximates the full 
NLO results quite well. 

Comparison of the NLO cross-sections with the cross-sections used in the 
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Figure 4: ii'-factors {K = (yNLo/<^Lo) foi' the four possible final states at (a): the Tevatron 
{m-g = 200GeV), and (b): the LHC {mq/m-g = 1.2). Parton densities: GRV94. 

experimental Tevatron analyses (LO, EHLQ parton densities, and a scale Q 
that equals the partonic eentre-of-mass energy) reveals that the NLO cor- 
rections raise the lower mass bounds for squarks and gluinos by -|-10GeV to 
+30GeV. For the LHC the shift in the mass values between LO and NLO 
cross-sections amounts to -l-lOGeV to +50GeV. 
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Figure 6: Normalized differential cross-sections for squark-gluino production at the LHC 
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Parton densities: 
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